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The Giant Radio Galaxy M 87 as a TeV 7-Ray Emitter 
observed with the HEGRA Cherenkov Telescopes 

Niels Getting, Martin Tluczykont and Gotz Heinzelmann 
for the HEGRA Collaboration 

Institut fiir Experimentalphysik, Universitdt Hamburg, Lumper Chaussee 149, 
D-22761 Hamburg, Germany 

Abstract 

For the first time an excess of photons above an energy threshold of 
730 GeV from the giant radio galaxy M87 has been measured with a signifi- 
cance of 4.1(7. The data have been taken during the years 1998 and 1999 with 
the HEGRA stereoscopic system of 5 imaging atmospheric Cherenkov telescopes. 
The excess corresponds to an integral flux of 3.3 % of the Crab flux. Making use of 
the imaging atmospheric Cherenkov technique, this is the first object of the AGN 
class observed in this energy range not belonging to the BL Lac type objects. 

1. Introduction 

Extragalactic TeV 7-ray emission has been observed with the imaging at- 
mospheric Cherenkov technique so far from AGN only of the BL Lac type, i. e. ob- 
jects ejecting matter in a relativistic jet oriented very close to the observer's line 
of sight. In BL Lacs, TeV photons are commonly believed to originate in the 
jets, most popularly due to inverse Compton scattering. The well studied ob- 
jects Mkn421 {z = 0.030), MknSOl {z = 0.034), lES 1959+650 {z = 0.047) and 
H 1426+428 {z = 0.129) belong to this type of TeV 7-ray emitters. However, 
other types of AGN, e. g. giant radio galaxies, also show jets, though, in contrast 
to BL Lac type objects, under large viewing angles. Amongst these the nearby 
radio galaxy M87 {z = 0.00436) - containing a supermassive black hole with 
Mbh ~ 2-3 X 10^ Mq [11] - has been speculated to be a powerful accelerator of 
cosmic rays, including the highest energy particles observed in the universe, see 
e.g. [10, 8]. The angle of the M87 jet axis to the line of sight was determined to 
be 30° - 35° [7]. 

The HEGRA collaboration has extensively observed M87 in 1998 and 
1999 with the stereoscopic system of 5 imaging atmospheric Cherenkov telescopes 
(I ACT system) [9] as one of the prime candidates for TeV 7-ray emission from the 
class of radio galaxies. In this paper the encouraging results of these observations 
are reported, applying a very sensitive analysis method (see also [3]). Astrophys- 
ical conclusions concerning the nature of the observed excess are discussed. 
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Fig. 1. Left: Number of events vs. squared angular distance to M87 as observed 
in 1998 and 1999 with the HEGRA lACT system (dots: ON-source events, dashed 
histogram: background). The significance of the excess amounts to 4.1a. Right: 
The center of gravity position of the HEGRA M 87 TeV 7-ray excess is marked by 
the cross indicating the statistical la errors (radio image adapted from [17]). 

2. Observations and results of analysis 

M87 was observed in 1998 and 1999 with the HEGRA lACT system for 
a total of 83.4 h [3] above a mean energy threshold of 730 GeV for a Crab-like 
spectrum [14]. Only data of good quality were considered for the analysis. 

All M87 observations were performed in the wobble mode allowing for si- 
multaneous estimation of the background ( "OFF" ) rate induced by charged cos- 
mic rays [1] . This analysis uses a ring segment as extended OFF-region reducing 
the statistical error on the number of background events [3]. The radius and width 
of the ring are set according to the position and size of the ON-source area in 
order to provide the same acceptance for ON- resp. OFF-source events. The event 
reconstruction (described elsewhere, e.g. [2]) makes use of algorithm #3 for the 
reconstruction of the air shower direction [13]. A "tight shape cut" (parameter 
mscw < 1.1) [14] is applied for an effective 7-hadron separation. The optimum 
angular cut was derived using 7-ray events from the Crab nebula on the basis of 
a nearly contemporaneous data set at similar zenith angles. 

Figure 1 (left panel) shows the event distribution for the ON-source and 
the OFF-source regions as a function of the squared angular distance of the shower 
direction to the source position. The statistical significance of the observed excess 
from M87 is 4.1 a, calculated using formula (17) from [16]. On the basis of 
the limited event statistics the excess is compatible with a point-like source for 
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the HEGRA lACT system at the position of M 87, although a sUghtly extended 
emission region cannot be excluded. 

The event distribution in the field of view was used to determine the center 
of gravity position (CoG) of the TeV 7-ray excess at aj2ooo.o = 12'"'30™54.4'' ± 
6.9L,±1.7:^^,, 5,,ooo.o = 12°24'17"±1.7L,±0.4;^, (see Figure 1, right panel). The 
accuracy of the CoG determination is limited by a systematic pointing error of 
about 25" [18]. Within the large statistical errors, the CoG is consistent with 
the M87 position. Currently, it is not possible to localize a candidate TeV 7-ray 
production site to particular inner radio structures of M87. 

The observed excess can be converted into an integral fiux of (3.3 ±0.8) % 
of the Crab nebula fiux. A conversion into absolute fiux units results in an M87 
7-ray fiux of N.,{E > 730 GeV) = (0.96 ± 0.23) x lO'^^ pj^Qt_ ^.^^-2 ^ spectral 
analysis of the data of the M87 data has been performed using the analysis 
technique described in [4]. The data can be well described with a power law 
dN/dE ~ with a = 2.9 ± 0.8,tat ± 0.08,,,^ 

3. Summary and Conclusions 

The radio galaxy M 87 has been observed with the HEGRA lACT system 
for a total of 83.4 h. For the first time a significant excess of 4.1 a has been 
detected at energies above a mean energy threshold of 730 GeV from a member 
of this class of objects using the imaging atmospheric Cherenkov technique. 

Due to the limited number of excess events detected so far it is difficult to 
draw a conclusion about the origin of the TeV 7-radiation. Assuming a spectral 
shape following a power law dN/dE oc E'"^'^ the integral photon fiux of (3.3 ± 
0.8) % of the Crab nebula fiux converts into an energy fiux of Ej{E > 730 GeV) = 
(4.3 ± 1.0) X lO-^^ergcm-^s^^ Given the distance to M87 of about 16Mpc, this 
corresponds to a 7-ray luminosity above 730 GeV of about lO^^ergs"^ under the 
assumption of isotropic emission. The integral fiux observed by HEGRA is not in 
confiict with an upper limit of 2.2 x 10~^^ phot. cm~^s~^ above 250 GeV reported 
by the VERITAS collaboration from data collected in 2000 and 2001 [15]. 

Several different possibilites for the origin of GeV/ TeV 7-radiation are 
conceivable. M 87 with its pc scale jet has recently been modeled within the Syn- 
chrotron Self Compton scenario as a BL Lac object seen at a large angle to its jet 
axis [5]. Note, that a recent Chandra monitoring of the optical knot HST-1 in the 
M87 jet (located only 0.8" away from the core) revealed a strong hint for a syn- 
chrotron origin of the observed X-ray emission [12]. The TeV 7-radiation of M87 
has also been modeled using the so-called Synchrotron Proton Blazar model [19]. 
In both models, the flux observed by HEGRA can be accommodated. The large 
scale (kpc) jets with several knots detected at radio to X-ray frequencies is also 
a possible 7-ray production site in M87. Moreover, 7-rays could be produced in 
the interstellar medium at larger distances from the center of M87. It should be 
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noted that M87 is also considered as a possible source of TeV 7-rays from the 
hypothetical neutralino annihilation process [6]. 

A weak signal at the centi-Crab level is at the sensitivity threshold for 
the HEGRA lACT system. Deep observations of M87 with next generation 
Cherenkov telescopes like H-E-S-S, MAGIC and VERITAS will provide a high 
sensitivity together with a low energy threshold. Due to the proximity to M87 
(16 Mpc compared to 110 Mpc for the closest TeV BL Lac Mkn421) these measure- 
ments will allow an accurate location and spectral analysis of the 7-ray emission 
site in M87 thus greatly advancing our understanding of its TeV 7-radiation. 
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Highlights from 6 years of TeV gamma-ray astrophysics 
with the HEGRA imaging Cherenkov telescopes 

Gotz Heinzelmann^ , for the HEGRA Collaboration^ 

(1) Universitdt Hamburg, Institut fiir Experimentalphysik, Luruper Chaussee 149, 
D-22761 Hamburg, Germany 

(2) http : / /www-hegra . desy . de/hegra/ 

Abstract 

The HEGRA (High Energy Gamma Ray Astronomy) experiment achieved out- 
standing results during the operation of the six Cherenkov telescopes (end of 1996 - end 
of 2002), pioneering with 5 telescopes the stereoscopic observation mode. Concerning 
the Galactic sources these include the detection of Cassiopeia A, being the only shell 
type supernova remnant at TeV energies seen up to now in the northern sky and re- 
cently the observation of a yet unidentified TeV gamma-ray source TeV J2032-I-4130 
in the Cygnus region. Also a scan of a large fraction of the Galactic plane has been 
achieved. Concerning the extragalactic sources most precise spectra at the highest en- 
ergies have been obtained from the well studied blazars Mkn 501 and Mkn 421 as well 
as from H1426-I-428 and lES 1959-1-650 established in the last two years only. Also ex- 
tensive multi-wavelength campaigns have been successfully performed. Recently strong 
evidence for the nearby giant radio galaxy M87 being a TeV 7-ray source has been 
obtained. Some of the highlight results are addressed below. 

1. Introduction and experimental method 

The imaging atmospheric Cherenkov telescopes (Table 1) of HEGRA (Ca- 
nary island La Palma, 28.75*^ N, 17.89" W, 2200m a.s.l.) started with a prototype 
telescope in 1992 and consisted in fall 1996 of a stand alone telescope (CTl) [1] 
pioneering observations during partial moon time and 5 telescopes (CT2-CT6) [2] 
introducing the very successful stereoscopic observation mode adopted by most of 
the next generation experiments. The stereoscopic observation of an air shower 
i.e. the simultaneous measurement of an air shower with several telescopes, under 
different viewing angles, allows to reconstruct unambiguously the shower direc- 
tion and impact point and also the shower height on an event by event basis. This 
leads to an improved angular and energy resolution, an improved gamma/hadron 
separation and suppression of background from night sky light and muons (due 
to the coincidence method). The stereoscopy also allows to record simultaneously 
events from well defined background regions giving the quoted errors a high cred- 
ibility and to perform sky searches in the field of view of the camera. The flux 
sensitivity achieved is a 10 a detection within 1 hour for a source with a flux of 
1 Crab. The telescopes CT2-CT6 ceased operation end of 2002. 
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Table 1. Properties of the HEGRA Cherenkov telescopes (end of 2002) 



Telescope 


Observ. 


Mirror 


PM-Camera 


FoV 


'^thres 


69 


6E/E 




mode 


area [m^] 


^pixels 


[deg] 


[GeV] 


[deg] 


[%] 


CTl 


alone 


5 ^10 


37 ^ 127 


3 


700-900 


< 0.2 


~25 


CT2-CT6 


stereo 


5 X 8.5 


5 X 271 


4.3 


500-600 


< 0.1 


10-20 



2. Galactic sources and Galactic accelerators of nuclei 

It is one of the main goals of the TeV 7-astronomy to find the sources of the 
Galactic cosmic rays commonly believed to be the shell type supernova remnants. 
In fact HEGRA detected [3] with Cassiopeia A the only shell type supernova 
remnant seen so far at TeV energies in the northern sky after a long observation 
time of 210 h (in 3 years) resulting in a fiux of ~ 3 % of the Crab nebula. Although 
not yet unequivocally proven it is very likely [4] that Cas A is one of the long 
sought hadronic accelerators. A scan of 1/4 of the Galactic plane (longitude 1= -2° 
to 85°) yielded upper limits in the range of 0.15-0.5 Crab fiux units for 63 SNR's, 
86 pulsars and 9 unidentified EGRET sources [5]. Recently HEGRA detected 
[6] with TeV J2032-I-4130 (Fig. 1) a new source up to now unidentified (i.e. no 
counterpart at radio, optical and X-ray energies) in a direction about 0.5 degrees 
north of Cygnus X-3. The source is detected meanwhile with a significance above 
7cr [6], exhibits a hard spectrum (index a ~ -1.9) and is possibly extended (on 
a 3o" level). Several source mechanisms are conceivable and it may turn out that 
this source plays an important role in the context of the Galactic accelerators. 



Table 2. Galactic and extragalactic sources seen by HEGRA. The fluxes are given in 



■2.59±0.03±0.05 



units of the Crab-flux [dJ^/dE = (2.79 ± 0.02 ± 0.5) • 10~^(E/lTeV)" 
photons m-2 s'^ TeV"!]. (C=CAT, H=HEGRA, V=VERITAS, T=7Tel. Array) 





Type 


distance 


si^ 


^nificance 


Flux 


Experim. 












(HEGRA) 


[Crab flux] 




Crab nebula 


Plerion 


1.6 kpc 


> 


10 


1.0 


many 


Cas A 


SNR (shell) 


3.4 kpc 




6 


0.03 


H 


TeV J2032+4130 


unknown 


unknown 




7 


0.03 


H 


Mkn 421 


BL Lac 


(z = 


0.030) 


> 


10 


0.04-7.40 


many 


Mkn 501 


BL Lac 


(z = 


0.034) 


> 


10 


0.33-6.00 


many 


lES 1959-^650 


BL Lac 


(z = 


0.047) 


> 


10 


0.05-2.20 


H,V,C,T 


H1426-h428 


BL Lac 


(z = 


0.129) 




7 


0.03-0.08 


H,V,C 


lES 2344-^514 


BL Lac 


(z = 


0.044) 




4 


0.03 


H,V 


M87 


radio gal. 


(z = 


0.0044) 




4 


0.03 


H 
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Fig. 1. Left: Skymap of the Cygnus region with the center of gravity and the 2a 
error circle for TeV J2032+4130 [6]. Also given are the 95% error ellipses of 
EGRET sources, the core of Cygnus 0B2 and the location of Cyg X-3. 
Right: Radio image of M87 at 90 cm showing the structure of the M87 halo. 
The center of gravity position of the TeV 7-ray excess from the HEGRA M87 
observations is marked by the cross (statistical \a errors) [11]. 

3. Extragalactic sources 

Right after the completion of 4 telescopes ready for the stereoscopic ob- 
servation Mkn 501 showed strong flares in 1997 during the whole observation 
period, allowing a measurement of an extended light curve and energy spectra of 
unprecedented precision in this field up to 16 TeV [7]. Remarkably, the shape of 
the energy spectrum essentially did not depend on the flux level in 1997. This is 
in contrast to Mkn 421 which exhibited also flaring states and showed examples 
for a hardening of the spectrum with rising flux level in the HEGRA data [8]. The 
Mkn 501 spectrum however became softer in the subsequent years. In the last 
two years two more blazars have now been well established, namely 1ES1959+650 
and H1426-I-428. 1ES1959+650, reported by the Seven Telescop Array with 3.9cr 
(1998), has been detected by HEGRA end of 2001 with 5.2 a and showed sub- 
sequently strong flares in May 2002 reaching a signiflcance level above 20cr [10] 
in the HEGRA data (detected also by VERITAS and CAT). H1426+428 is the 
most distant TeV 7-source up to now and has been detected by all three northern 
instruments VERITAS, HEGRA and CAT (Table 2). 

The precise HEGRA energy spectrum of Mkn 501 has been used by sev- 
eral authors to infer the density of the extragalactic background light (EBL) in 
the optical to infrared range due to the absorption process ^Tev + 1ebl e~^e~. 
In fact, the spectra of Mkn 501 (and Mkn 421) in the TeV range yield the best 
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limits for the EBL density in the mid-infrared range to date, which is of great 
cosmological interest. With more precise data also from other sources such as 
H1426+428, which is a factor four more distant (and with a better understanding 
of the sources) even more stringent results will be possible. The present HEGRA 
data for H1426+428 already indicate a modulation of the measured energy spec- 
trum [10] due to the imprint of the absorption process mentioned above. To 
understand the underlying acceleration processes in the jets of the blazars several 
multi-wavelength campaigns have been performed e. g. together with the X-ray 
satellite RXTE especially for Mkn 421 and Mkn 501. Correlated flux variations 
have been observed which indicate that an inverse Compton origin is in general 
able to describe the TeV data though other processes can not be excluded. 

Analyzing the data from observations of further 52 objects from the AGN 
class with a total observation time of 658 h yielded signals on the 4 a level from 
the blazar 1ES2344-I-514 [11] (reported earlier by Whipple) and very recently from 
the radio galaxy M87 (Fig. 1) [12] which still has to be confirmed. M87 at the 
center of the Virgo cluster is of specific interest since it is the first TeV AGN being 
observed using the imaging technique and not belonging to the blazar class, thus 
opening the field of TeV 7-astronomy to a new class of sources. This very nearby 
object may also play an important role for the acceleration of cosmic rays up to 
the highest energies. Further results will be presented at this conference. 

4. Conclusion 

The contributions of HEGRA to the TeV 7-astronomy on the way to a 
mature field have been extremely important by pioneering new experimental tech- 
niques and related analysis methods, by obtaining high statistics and precise data 
for light curves and energy spectra and by detecting new sources such as Gas A, 
TeV J2032+4130, lES 1959+650 and M87, thus promising a great future to the 
next generation Gherenkov telescopes. 
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Search for TeV Gamma- Rays from the Andromeda Galaxy 
and for Supersymmetric Dark Matter in the Core of M31 

W. Hofmann, D. Horns and H. Lampeitl, for the HEGRA collaboration 
Max-Planck-Institut fiir Kernphysik, D 69029 Heidelberg, P.O. Box 103989 

Abstract 

Using the HEGRA system of imaging atmospheric Cherenkov telescopes, 
the Andromeda galaxy (M31) was surveyed for TeV gamma ray emission. Given 
the large field of view of the HEGRA telescopes, three pointings were sufficient 
to cover all of M31, including also M32 and NGC205. No indications for point 
sources of TeV gamma rays were found. Upper limits are given at a level of a few 
percent of the Crab flux. A speciflc search for monoenergetic gamma-ray lines 
from annihilation of supersymmetric dark matter particles accumulating near the 
center of M31 resulted in flux limits in the 10~^^ cm~^s~^ range, well above the 
predicted MSSM flux levels except for models with pronounced dark-matter spikes 
or strongly enhanced annihilation rates. 

1. Introduction 

M31, as the nearest large galaxy at a distance of 770 kpc, most likely 
with a black hole of 3 ■ 10'' M© at its core [7], is an obvious target for TeV 
observations. While conventional sources such as the Crab Nebula would not be 
visible over the distance to M31, given the current sensitivity of the instruments, 
the energy output of some of the X-ray sources detected in M31 is at a level which, 
if continued to higher energy, might be detectable. A particular mechanism for 
gamma-ray emission is the annihilation of supersymmetric dark matter particles 
accumulating at the core of M31; rotation curves suggest that M31 contains a 
signiflcant amount of dark matter. 

2. The HEGRA telescope system and the M31 data set 

M31 was observed in August, September and November 2001, using all 
five telescopes of the HEGRA system of Cherenkov telescopes on the Canarian 
Island of La Palma. Each of the telescopes has a 8.5 m^ mirror and is equipped 
with a 271-pixel photomultiplier camera with a 4.3° field of view. The system has 
an energy threshold around 500 GeV and a 0.1° angular resolution. The detection 
rate for gamma rays is rather uniform within the central 2° of the field of view, 
and drops to 63% of its peak value at 1.8° from the optical axis. In order to cover 
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all of M31, observation time was distributed between three tracking positions, one 
centered on the core of M31 at RA Oh 42' 44" DEC 41° 16' 9.12", one displaced 
by 0.56° to the SW at RA Oh 40' 30" DEC 40° 39' 0" and one displaced to 
the NE at RA Oh 44' 43" DEC 41° 53' 0". After data cleaning, 20.1 h of good 
observation time were selected, most of which were taken at Zenith angles below 
25°. For calibration and reference, 9.7 h of Crab Nebula data taken in October 
and November were selected. 

In a first analysis step, the entire field was searched for indications of 
gamma-ray point sources. The data analysis and selection of gamma-ray candi- 
dates follows to a large extent the procedures developed for earlier surveys [1,2]. 
For each grid point, gamma ray candidates with reconstructed directions within 
0.143°, consistent within the angular resolution, were counted. Since no specific 
off-source data were taken, the expected background for each grid point is deter- 
mined from the data itself, using appropriate background regions in equivalent 
locations of the field of view of the cameras. 

No obvious excess of signal events is visible at any grid point. The sig- 
nificance for each grid point is then calculated based on [8]. The distribution of 
significances for all grid points approximates a Gaussian distribution with unit 
width, as expected in the absence of a genuine signal; the point with the highest 
significance of 3.5 a is well compatible with the expected tail of the Gaussian 
distribution, and should be interpreted as an upward fluctuation. 

We conclude that there is no statistically significant indication of a TeV 
gamma ray point source in M31, and extract upper limits for the source flux. The 
99% confidence level on the number of excess events is calculated for each grid 
point, and is converted into a flux using the measured gamma-ray rate from the 
Crab Nebula. The resulting flux limits are shown in Fig. 1 (left), and range from 
3.3% of the Crab flux for the center of M31 to about 30% at the periphery. For 
two other objects in the survey range, M32 and NGC205, limits of 4.4% and 2.8% 
of the Crab flux are derived, respectively. More details can be found in [3]. 

3. Search for supersymmetric dark matter in M31 

Beyond the general survey, the good energy resolution of the HEGRA 
telescope system allows a dedicated search for supersymmetric dark matter in 
M31, looking for the line emission from the self-annihilation of the lightest stable 
particle (LSP). We have searched for line emission in the energy spectrum from the 
central part of M31 (r < 1.4 kpc), and compare the results with predictions of the 
flux of gamma-rays from the annihilation of neutralinos (xo) in the framework of 
the minimal supersymmetric standard model (MSSM) using a spatial distribution 
of the radial density of neutralinos in M31 in concordance with the measured 
rotation curve [6]. 

In order to reduce the background from isotropic cosmic ray events, in the 
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Fig. 1. Left: Upper limits on the flux from TeV point sources in M31, expressed 
in units of the Crab flux. For reference, the locations of the center of M31, the 
10 kpc dust ring and the two companion galaxies M32 and NGC205 are indicated. 
Tracking positions are indicated by crosses. Right: Exclusion region (grey shaded 
region) for the the gamma-ray line flux from the center of M31. The two sets of 
curves labeled 77 and Z^^ indicate the upper range of model predictions for the 
minimal supersymmetric models calculated in [4]. 



event selection a tight cut on the direction has been chosen for the central part of 
M31. Only events from within a cone of radius 0.105° are accepted. At a distance 
of 770 kpc this corresponds to the inner 1.4 kpc of M31. The sensitivity for line 
emission depends on the energy resolution of the detector. For this analysis, an 
improved energy reconstruction algorithm has been applied. The relative energy 
resolution (AE/E) reaches 10% for a wide energy band from the threshold of 
500 GeV up to 10 TeV. The search bin in energy is 12 % wide, concentrating on 
the central section of a monoenergetic line in order to achieve best signal to noise 
ratios. The energy bins are equally spaced on a logarithmic scale such that the 
bin centers are separated by 0.025 on a decadic logarithmic scale, making the bins 
correlated. 

The background expectation of the measurement is determined using seven 
independent OFF-regions with similar acceptance to the ON-region. For the cal- 
culation of upper limits on the number of excess events N^^'^\ the 99 % c.l. upper 
limits were caculated and an upper limit on the rate from a 7-ray line was derived. 
The corresponding flux limits were calculated using collection areas A^ff^E, 9) de- 
rived from Monte Carlo simulations applying the same reconstruction methods 
and event selection as for the data analysis. The resulting exclusion region is 
indicated in Fig. 1 (right) as the grey shaded region. 

A prediction for the line flux emission was derived using scans of the 7- 
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parameter space of MSSM performed by [4] . The expected flux is given by Eq. (6) 
of [6] or the equivalent Eq. (13) of [4] (scaled to the distance to M31): 

= (3.18 ■ 10^'' photons cm-2 s"!) ( 1'^''^^^ 

\10~^° cm'^ s 

Apart from trivial factors, it depends on the velocity-weighted annihilation cross 
section (erf) and on the line-of-sight and solid-angle integrated squared LSP mass 
density p 

S = J J dsdn 

written as Sig in units of 10^^ GeV^ cm~^. The envelope of the allowed (av) 
vs. space was used as given in Figs. 1 and 2 of [4]. The unitless parameter 
Si9 ^ 1.7 corresponds to Sig = 3 of [6], scaled to the inner 1.4 kpc on the basis 
of the p(r) distributions of Fig. 3(b) of [6] and using a distance of 770 kpc. 
This value of Sig represents the upper range of the different models for the dark 
matter halo of M31 discussed in this reference. The predicted line flux is given in 
Fig. 2 (right). We have indicated the two individual final states with line emission 
^ 77 ^iid ~^ ^^l) ill the figure. The model prediction given here is 
based upon smoothly distributed dark matter. Within this model, no signal from 
M31 with the current sensitivity of Cherenkov detectors is to be expected. Very 
favorable conditions as for example dumpiness of dark matter would lead to a 
considerable increase of the annihilation rate [5,9]. It is conceivable that the flux 
level of some models would become detectable by tuning the distribution of the 
dark matter halo and by invoking mechanisms to allow for larger cross sections. 
However, in the absence of a signal, no exotic speculations are required. 
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Studies of the Crab Nebula based upon 400 hours of Obser- 
vations with the HEGRA System of Cherenkov Telescopes 

D. Horns {dieter.horns@mpi-hd.mpg.de), for the HEGRA collaboration 
MPI f. Kernphysik, Postfach 103980, D-69 117 Heidelberg, Germany 

Abstract 

The Crab nebula has been observed extensively with the HEGRA system 
of imaging air Cherenkov telescopes. Roughly 400 hours of prime observation 
time at zenith angles between 5 and 65 degrees have been dedicated to the stan- 
dard candle in TeV astronomy. Based upon the data set gathered during the 
5 years of operation of the HEGRA telescopes, the energy spectrum has been 
reconstructed. The energy spectrum extends from 500 GeV beyond 20 TeV and 
allows to constrain the position of a possible high-energy cut-off. 

1. Introduction 

The Crab nebula has been a prime target of observations since the begin- 
ning of operation of the instrument primarily for calibration and monitoring of the 
instrument's 7-efficiency in the initial stage of operation and during the routine 
operation. The observational results include the energy spectrum [3], upper lim- 
its on the angular size of the emitting region [4], and constraints on the fraction 
of pulsed radiation [5]. After the end of the live-time of the instrument (it has 
been dismantled in October 2002), an analysis of the entire dataset taken from 
the Crab nebula with increased statistics and improved reconstruction technique 
with respect to previous publications [2] (the most recent published spectrum 
from the Crab nebula is based upon ^ 150 hrs of observations) is justified and 
first preliminary results on the energy spectrum are presented here. 

2. Observations and Results 

The observations had been carried out between September 1997 and Septem- 
ber 2002. During it's life-time as a 4 (later 5) telescope system, the setup of the 
individual telescopes (including cameras and electronics) has not been modified. 
However, the degradation of the mirrors and the ageing of the photomultipliers 
have been compensated by an increase of the high voltage on an annual basis. 
The performance of the instrument has been closely monitored by various checks 
on the acceptance, absolute calibration with muon-rings, cut efficiencies, point- 
ing calibration, and careful scrutinizing of the photomultiplier characteristics [7]. 
The benefit of the thorough calibration and understanding of the instrument's 
performance is the reliable data-taking that in turn allows to combine data taken 
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Table 1. Summary of observations, split up into 4 bins of altitude. 



Season 


alt= 65 - 


- 85° 50 - 65° 


40 - 50° 


15 - 40° E 


E 


year 


[ksec 


[ksec] 


[ksec 


[ksec 


[ksec] 


[hrs] 



1997/98 162.04 96.83 54.48 49.67 363.02 100.84 

1998/99 217.84 60.82 92.15 115.72 486.52 135.15 

1999/00 108.89 9.10 0.66 0.15 118.80 33.00 

2000/01 129.09 39.59 17.02 2.96 188.66 52.40 

2001/02 125.36 52.69 32.90 13.62 224.58 62.38 



over many years with a minimum of systematics introduced by changes of the 
detector. 

The simulations of the detector have been performed using two indepen- 
dent Monte-Carlo type air-shower and detector simulations [8]. Both simulations 
agree within 10 % on the relevant parameters (cut efficiency, rate of cosmic ray 
events, angular resolution, energy resolution). All detector simulations have been 
carried out for individual periods (one moon cycle) taking the degrading photo- 
multiplier and optical gain into account to model the response of the detector 
and particularly its change with time. Additionally, whenever required, the sim- 
ulations have been done with different telescope setups (3, 4, and 5-telescope 
setups) in case of failures of individual telescopes. The observations are taken in 
wobble-mode, where the telescopes point with an offset of 0.5° to the direction of 
the source, alternating the sign of the displacement from run to run. A summary 
of the observational time is given in Table ^ The listed observational times are 
selected for good weather conditions and instrument performance. Overall, less 
than 10 % of the runs have been rejected by requiring the event rate to deviate 
by less than 25 % from the expected value (for a given observational period and 
zenith angle). 

The analysis technique is based upon the method first introduced for 
HEGRA data in [2] with the modification of applying a slightly tighter cut on 
the main 7/hadron separation quantity mean scaled width< 1.1. The mean scaled 
width is related to the width of the images scaled to the expectation for a 7-ray 
induced event. 

For individual events at a given zenith angle, an energy estimate is calcu- 
lated with a relative resolution of AE/E < 12 % [6]. Using Monte-Carlo calcu- 
lated tables of collection function of reconstructed energy for 5 discrete 
zenith angles a collection area for each event is calculated using interpolation in 
zenith angle and energy. Finally, the sum of the inverse of the collection areas 
is calculated for a given energy bin and using the dead-time corrected on-time, a 
differential flux is calculated. The background is calculated by averaging over 5 
regions in the camera with the same distance to the camera center and diameter 
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as the source region. 

We have omitted to include the observations below an altitude of 40°. Fur- 
ther results on large zenith angle observations will be presented at the conference. 
The preliminary result is presented in Fig. El The deep exposure is sufficient to 
collect ~ 10 000 photons used in the analysis covering 2 decades in energy. The 
signal in the highest energy bin centered on 39 TeV is beyond 5 a. Inclusion 
of large zenith angle observations will overcome the difficulties of saturation and 
boost the statistics of photons up to even higher energies. 

The data is well described by a pure power-law fit dN^/dE = (2.91 ±0.08) ■ 
10-11 (E/TeV)-2-59±o-02 ph/(cm2 s TeV), xledi.d.o.f .) = 1.6(11). A functional form 
with a curvature term describes the data slightly better: dN^/dE = (2.72 ± 
0.08) ■ 10-" ■ (E/TeV)-2-43±o.o5-(o.i4±o.o4).iogio{i?/TeV-)^ xled{d.o.f.) = 0.77(10). A 

similar xledi.d-o.f) = 1(7) is reached for a power-law with an exponential cut-off: 
dN^/dE = (2.73 ± 0.08) ■ lO-^^ ■ (E/TeV)-2-43±o.o5 . exp(-E/(53 ± 38)TeV). The 
very slight curvature changes the photon-index from 2.43 at 1 TeV to 2.65 at 
40 TeV. 
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Fig. 1. The differential photon energy spectrum for the Crab nebula. The preliminary 
result is shown with the statistical errors only. The systematic errors are dominant 
below 1 TeV and above 20 TeV and are not included. Additionally, the measurement 
is subject to an overall uncertainty of the absolute energy scale estimated to be less 
than 15 %. 
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3. Discussion 

The observations presented here cover for the first time a wide range in 
energies (two decades) with a single instrument. The advantage of such an ob- 
servation is clearly that relative calibration uncertainties between different mea- 
surements are absent. However, there are systematic uncertainties that need to 
be addressed before conclusions on the physics are derived. The statistical errors 
of the measurement are very small (< 5 % for energies between 1 and 4 TeV) and 
therefore the systematic uncertainties are dominating. 

The strong energy dependence of the response function of the telescopes 
in the threshold region is introducing uncertainties in the region below 1 TeV. 
Different reconstruction methods and event selections have been carried out and 
the variation of the response at lower energies has been studied. The conserva- 
tive estimate of the systematic uncertainty below 1 TeV is a 50 % uncertainty 
region. This estimate is largely based upon the remaining variations in the en- 
ergy spectrum below 1 TeV in different zenith angle bands and over the years of 
observation. 

For the high energy end of the spectrum, saturation of the electronics 
used for digitization of the pulses is crucial. The dynamical range of the FADC 
electronics is limited to one order of magnitude before the individual time slices 
of the FADC saturate. However, using the overall pulse shape makes a larger 
dynamical range accessible. This introduces possible systematic effects by either 
over- or under-compensating. From the measured image amplitudes of cosmic 
ray induced events we expect these effects to be small (< 10 %). In a previous 
analysis [2] of large zenith angle observations where the image amplitude even 
for events with energies beyond 10 TeV are not saturating, no strong systematic 
effect was visible. 

Further studies on the pulsed emission from the Crab pulsar and morphol- 
ogy of the 7-ray emitting region at different energies are on the way. 
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High Energy Emission from H1426+428 and Absorption on 
the Extragalactic Background Light 

D. Horns, F.A. Aharonian, L. Costamante, for the HEGRA Coll. 
MPI f. Kernphysik, Postfach 103980, D-69117 Heidelberg, Germany 

Abstract 

The Blazar H1426+428 {z = 0.129) has been observed with the HEGRA 
stereoscopic system of imaging air Cherenkov telescopes for 270 hrs in the years 
1999, 2000, and 2002. The object is detected with a significance of 7.5 a. The 
energy spectrum measured by HEGRA extends up to ~ 8 TeV with indications 
for a change in the spectral slope at energies above 2 TeV. The integral flux above 
1 TeV in 1999 and 2000 amounts to ^ 8 % of the flux of the Crab Nebula. In 
2002, the flux drops to a level of 3 % of the Crab Nebula. The extragalactic back- 
ground light (EBL) at wavelengths between 1 and 10 microns causes substantial 
(r = 1 ... 9) absorption of photons in the measured energy spectrum. The in- 
trinsic spectra after correction for absorption are discussed. Nearly simultaneous 
observations with the X-ray observatory RXTE in 2000 and 2002 indicate flux 
and spectral variability in the X-ray band between 2 and 20 keV. 

1. Introduction 

The ongoing exploration of Blazar type objects over a broad region of the 
electromagnetic spectrum has revealed surprising properties that are of relevance 
for the study of injection, acceleration, and radiative cooling of particles in rel- 
ativistic jets. However, a deep understanding of these objects based upon the 
observations at TeV energies suffers from the ambiguity of the observed spectrum 
with respect to absorption effects on the extragalactic background light (EBL). 
Currently, the object H1426+428 at a red shift oi z = 0.129 is the most distant 
emitter of photons up to energies of ~ 8 TeV detected. Unfortunately, the ob- 
servational uncertainty on the measured energy density of the EBL is quite large 
(up to a factor of 10 at A ~ 1 /im, see also Fig. lb). Given the consequently 
large uncertainty on t{E), a reasonably accurate (better to within a factor of 10) 
correction for the effect of absorption is not possible. Here, another approach is 
followed: Using widely different descriptions of the EBL, absorption corrections 
are applied to the observed energy spectrum to infer the intrinsic source spec- 
trum. The intrinsic source spectrum is then checked for consistency with our 
current understanding of Blazar physics. 
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Table 1. Table of X-ray observations. The Galactic column density has been 
kept fixed at uh = 1-38 • 10^'^ cm~^. The fit range for a power law 
dN/dE = iVo(^/keV)-T starts at 4 keV. 



Pointing 


Exposure 


No 


7 


xVd.o.f.{d.o.f.) 


[year] 


[ksec] 


[keV-^s-^cm-^] 






2000 


12.7 


(9.1 ±9,7) • 10~3 


1.96 ± 0.04 


2.82(7) 


2001 


4 


(30.1 ± 1.6) • 10-3 


1.81 ±0.03 


0.75(23) 


2002 


110 


(9.5 ±0.3) • 10-3 


1.82 ±0.02 


3.76(27) 



2. Observations and Results 

2.1. X-ray observations 

The instruments on-board the Rossi X-ray Timing Explorer (RXTE) have 
been used to measure the X-ray flux from H1426±428 during the years 2000, 
2001, and 2002. The overall flux level is quite different between the observations, 
varying by a factor of 4 between 2000 and 2001 (see Table El for a summary of 
the pointed observations). During all three pointings a hard spectrum described 
by a power law with dN/dE oc E~'^ with 7 < 2 is observed. This indicates that 
the object is showing only weak spectral variations, remaining in a hard state for 
different flux levels. Interpreting the corresponding spectral energy distribution 
in the framework of a combined synchrotron and Compton dominated emission 
model, this would imply that the position of the synchrotron peak is at or beyond 
10 keV, and possibly similar to previous observations [3] located in the hard X-ray 
band. In 2000, the spectrum seems to soften above 8 keV. In 2002 the spectrum 
seems to harden beyond 15 keV. A pure power-law does not satisfactorily describe 
the data in these cases. 

The all-sky monitor (ASM) on-board the RXTE satellite offers a continu- 
ous coverage of the object. Averaging the daily measurements of the flux for the 
HEGRA observational nights during the 1999, 2000 data set results in an ASM 
rate of 0.33 ± 0.08 counts/sec whereas in 2002 the ASM count rate drops by a 
factor of 1.7 to 0.19±0.04 counts/sec. In the context of a leptonic emission model, 
a correlation of the X-ray flux and the TeV flux is expected and is confirmed by 
the HEGRA data (see next section). 

2.2. TeV observations 

The observations carried out with the HEGRA system of 5 imaging air 

Cherenkov telescopes detected the object in 1999/2000 during a 40 hrs observation 
at a flux level of 8 % of the Crab Nebula [1]. In order to verify the observed 
energy spectrum with better statistics, an extended observational campaign has 
been carried out in 2002, accumulating 217 hrs of good data [2]. However, the 
source showed a lower flux of merely 3 % of the Crab Nebula. The spectral 
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Fig. 1. a) The differential energy spec- 
trum of H1426+428 as measured by 
HEGRA in the two different flux states 
seen in 1999, 2000, and 2002. For com- 
parison, the previously published results 
are included, b) A selection of measure- 
ments and models for the EBL from UV 
to far infrared (taken from [2]). 



analyses presented here are based upon an improved analysis technique tailored 
for weak sources. The event-selection criteria have been adopted to weak sources 
by applying tighter cuts on the arrival direction and on mean scaled width in 
order to increase the signal-to-noise ratio. An improved energy reconstruction 
algorithm has been used which benefits from the stereoscopical reconstruction 
of the position of the shower maximum and achieves a relative energy resolution 
(AE/E) of 15 % at the threshold and of 10 % at higher energies. The new method 
has successfully been used for data taken on the Crab Nebula and on the original 
40 hrs data set of H1426+428 taken in 1999 and 2000. The resulting data points 
are shown in Fig. ^ Generally, good agreement between the previously published 
energy spectrum and the reanalysed data of 1999/2000 is seen. The 2002 data 
confirm a hardening in the energy spectrum above 1 TeV. Combining the excess 
above 2 TeV of all the data results in a significance of 6 a. 

3. Conclusions 

Based upon the fact that similar flux levels of H1426+428 have been ob- 
served by the CAT, Whipple, and HEGRA group and that the ASM observations 
confirm a similar level for these observations at X-ray energies, it seems feasible 
to combine the different data sets to benefit from a wider energy band covered. 
Especially at energies below 1 TeV, a steep rise of the energy spectrum because of 
absorption effects is expected and confirmed by the different groups [4]. The re- 
sulting energy spectra are shown in Fig. 2(a-c) together with absorption corrected 
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Fig. 2. The differential energy spectrum combined from the measurements of the 
CAT, Whipple, and HEGRA (1999/2000) groups multiphed by E"^ (filled symbols). 
Prom left to right (a-c): The observed and corrected spectra for three different mod- 
els of the extragalactic background light (PrimOl, Aha02, MSOl, see also Pig. lb). 
The open symbols indicate the corrected spectrum. The solid lines are power law 
and broken power-law fits to the open symbols. The dashed lines are the same 
power-laws after multiplying the absorption term exp(— r(£')). 



spectra using the SEDs of the EBL shown in Fig. lb. Note, the models Aha02 
and MSOl result in an unreasonable rise of the source spectrum {dN/dE oc E) 
either at the low energy part (Aha02, Fig. 2b) of the spectrum or the high energy 
part of the spectrum (MSOl, Fig. 2c). Consulting Fig. lb it becomes clear that 
the Aha02 description of the EBL data (especially the high flux between 1-2 /im) 
results in a steep rise of r for E > 100 GeV. The high value of the SED in the 
near-infrared as given by MSOl causes a quickly increasing optical depth above 
2 TeV. 

In summary, the TeV observations of extragalactic objects offer an indirect 
approach to constrain the EBL. The HEGRA observations of H1426-I-428 with a 
clear detection of 6 a above 2 TeV are important to constrain the shape of the 
SED of the EBL in the near infra-red above 2 fim. Combined data from CAT, 
Whipple, and HEGRA constrain the EBL at wavelengths below 2 /xm. 
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An new method to determine the arrival direction of indi- 
vidual air showers with a single Air Cherenkov Telescope 

Daniel Kranich/'^ and Luisa Sabrina Stark, ^ for the HEGRA Collaboration 

(1) Max-Planck-Institut fiir Physik, Miinchen, Germany 

(2) Swiss Federal Institut of Technology Zurich, Zurich, Switzerland 

Abstract 

We present a new method to reconstruct the arrival direction of individual 
air showers. The method is based in part on the arrival direction reconstruction 
method of Lessard et al. [1] from the Whipple collaboration, but yields a sig- 
nificant 30% improvement in the obtained angular resolution. The method was 
successfully tested on Monte Carlo Simulations and real data of Mkn 421, Mkn 
501 and 1ES1959 from the HEGRA Cherenkov Telescope CTl. Based on the 
same data an angular resolution of 0.1° could be derived. 

1. Introduction 

The determination of the arrival direction of 7-ray photons is essential for 
the investigation of extended sources and the search for sources with unknown 
or inaccurate positions. In the case of an observation with a few Cherenkov 
telescopes the images of at least 2 (sometimes 3) detectors are needed in order 
to determine the arrival direction without ambiguity. However, even for a single 
telescope it is possible to circumvent this problem. Since all 7-rays have a similar 
orientation (limited by the source's extension) one can use the images of several 
shower events to estimate the arrival direction of an individual 7-ray photon. 

2. The Method 

The method to derive 2-dimensional source maps for a single telescope is 
based on the procedure described in [1] and which can be summarized as follows: 
The arrival direction for a given shower event is estimated as the point on the 
major axis of the shower image located at a distance DISP to the shower centroid. 
The DISP parameter is a function of the elongation of the shower image and is 
defined as DISP = ^ (1 - WIDTH/ LENGTH). The scaling parameter ^ has to 
be determined from real observations or MC data. The final source map is build 
up from the arrival directions of all shower events. In the case of the CTl data, 
the DISP parameter was slightly modified to include the LEAKAGE parameter: 
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DISP = ^ 1 



LENGTH • (1 + r; • LEAKAGE) 



WIDTH 



) 



(1) 



The LEAKAGE parameter is defined as the ratio between the hght content in the 
camera edge pixels and the total light content of the shower event. The inclusion 
of this parameter corrects for truncation effects of the LENGTH parameter in the 
small CTl camera ( 3° diameter). The free variables in Equ. ^ ^ and r] were 
determined from Monte Carlo simulations as ^ = 1.3 and rj = 6.6. 

As a modification of this method we determine a set of possible arrival 
directions for each shower event (in the following called arrival distribution): 
Taking into account the value and error of both the DISP parameter and the 
orientation of the major axis of each shower image one can calculate the most 
probable intersection point for any (randomly chosen) triple of shower images by 
means of a fit (see Fig. El)- The arrival distribution for a given shower event 
is then defined by the subset of all those intersection points, where this event has 
been part of the triple*. Poorly defined intersection points with a x^/dof. > 6/3 
were rejected. Each derived arrival distribution was then normalized to unity. The 
final ON source map is given by the superposition of arrival distributions from 
all individual shower events and the final excess arrival distribution is derived by 
subtracting the normalized OFF from the ON source map. 

The basic idea behind this approach is that all intersection points of pure 7 shower 
triples should point towards the same sky region whereas a more or less isotropic 
distribution is expected whenever hadronic shower images are involved. 

3. Results and Conclusions 

The excess arrival distribution for the Mkn 421 CTl data sample is shown 
in Fig. ED A fit of a two dimensional Gaussian onto the data has been used to 
derive estimates for the angular resolution {a^A and a dec) and the reconstructed 
source position (i?y4rec and DEGrec)- The results from this fit on several CTl 
data samples are shown in Table IT1 for both, the original method of Lessard et al. 
[1] and the modified method as presented here. As can be seen from the table, 
the new method yields a significant 30% improvement in the angular resolution 
for almost all sources. The reconstructed direction of the objects coincides with 
the known position within the accuracy of the shaft encoder steps (0.02°). It 
is also evident from Table El that the angular resolution strongly decreases for 
off-axis sources (i.e. where the telescope was pointing aside the source). This 
effect seems not to be related to the method itself since data which was artifi- 



*Due to the large computational overhead of this method the number of intersection points 
per shower event have been limited to 500. 



27 




gives the la error of the intersectionFig. 2. The excess arrival distribution for 
point. The 2 dimensional distribution the Mkn 421 CTl data sample, 
denotes the final arrival distribution for 
the black shower event. 

cially shifted towards an off-axis position^ didn't result in a decreased angular 
resolution. One possible reason could be that for off-axis observations the small 
CTl camera (about 3° diameter) leads to a larger number of truncated images, 
which then lower the effectiveness of both arrival direction estimation methods. 
This is currently under investigation. 

The decrease in angular resolution between the 1997 (Mkn 501) and the 2001/2002 
(Mkn 421, 1ES1959) data is addressed to the CTl mirror upgrade at the end of 
1997. The new CTl mirror was increased from 5 m^ to 10 m^ and shows some 
increased aberration effects due to the larger mirror diameter. 
Even though the derived angular resolution for 1ES1959 is similar for both meth- 
ods, the distribution of the excess arrival direction as obtained with the new 
method is superior (see Fig. E]) and yields a much clearer excess peak. 

In summary the presented method to determine the arrival direction of 
individual photons for a single Cherenkov telescope gives a 30% improvement 
compared to the original method. The derived angular resolution of ~ 0.1° is 
similar to the result one would obtain from a system of Cherenkov telescopes 

^Here the image parameters ALPHA and DIST were recalculated relative to a shifted camera 
center while the WIDTH and LENGTH parameters were kept. 
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Fig. 3. The excess arrival distribution for the 1ES1959 data sample of 2002 as deter- 
mined with the new (left) and original (right) method. 



Table 1. Results of the new (this paper) and original [1] arrival direction reconstruc- 
tion method as derived from different CTl data samples, ctra and udec denotes the 
angular resolution in RA and DEC direction, RA^qc and DECxec the correspond- 
ing reconstructed source position. Ara and Adec denote the difference between 
the real and the reconstructed source position {Ara ■= -Resource — -R^rec and 
Adec '■= DECsouice — DECxec)- The term 'oa' denotes data where the telescope 
was not directed towards the source but at an angular distance of 0.3°. 



data sample 




new 


method 






original method 






O'DEC 


Ara 


Adec 




O'DEC 


Ara 


Adec 


Mkn 421, 2001 


0.075 


0.073 


0.005 


0.014 


0.100 


0.099 


0.004 


0.015 


Mkn 501, 1997 


0.058 


0.056 


-0.023 


-0.035 


0.078 


0.077 


-0.023 


-0.035 


Mkn 501, oa, 1997 


0.082 


0.068 


0.017 


-0.010 


0.109 


0.101 


0.008 


-0.030 


1ES1959, 2002 


0.077 


0.085 


0.035 


-0.001 


0.082 


0.088 


0.052 


0.003 


MC 


0.052 


0.052 


-0.010 


0.022 


0.073 


0.072 


-0.010 


0.023 


MC oa 


0.081 


0.081 


-0.040 


0.019 


0.119 


0.123 


-0.039 


0.022 



[2]. This improvement should also allow for an increased cut sensitivity, which, 
however, has not been investigated yet. 
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Abstract 

Between 1997 and the end of the detector operation in fall 2002, about 5500 
hours of observational data were recorded with the HEGRA system of Imaging 
Atmospheric Cherenkov Telescopes (lACT). Besides dedicated scan observations 
of extended sky regions, a considerable fraction of the sky has been looked at as 
a side effect during the observations of selected source candidates. Altogether, 
more than 3% of the total sky has been observed with the HEGRA lACT system. 
We report on a search for possible new TeV sources within this entire data set. 

1. Introduction 

A large fraction of the observation time of the HEGRA lACT system [3] 
was dedicated to individual objects such as known supernova remnants and active 
galactic nuclei, which were known from other wavebands. Given the system's large 
field of view (FoV) of 4.3° with a homogeneous 7-ray acceptance over more than 2° 
in diameter, scans of extended sky regions were also possible (e.g. [2]). Summing 
up scans as well as observations of individual objects, more than 3% of the entire 
sky has been observed with the HEGRA telescope system. 

In the case of observations which were targeted at a point source, a large 
fraction of the FoV is not affected by the possible emission from the target object, 
and is usually only used to derive background estimates for the source candidate. 
It is however obviously possible, in analogy to scanning observations, to regard 
any position in the FoV as possible source candidate, and derive background 
estimates from other parts of the FoV. 

2. Analysis challenges 

The main analysis challenge is a reliable background estimate for an ar- 
bitrary position in the FoV. Gamma-ray event candidates are identified against 
the much larger hadronically induced background solely by the shapes of the 
event images. The image cut typically rejects 92% of the hadrons, the remain- 
ing hadrons generate background in the 7-ray shape cut domain. Depending 
on the total observation time and the angular extension of the search area, the 
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Fig. 1. Acceptance change in the field of view due to stars. The figures show the 
acceptance relative to a ring segment around the star positions, as a function of 
the blue star magnitude. The data points were derived on average for all stars in 
magnitude intervals. At the left panel, 7-ray like events were used, on the right 
side hadron-like events. The different symbols indicate different zenith angle bands. 
The lines are parametrisations which are used to account for the acceptance change. 



systematic accuracy of the determination of the absolute background level at the 
source candidate location must be 1% for very deep observations. 

In HEGRA system data analysis, the standard approach to determine a 
background estimate uses only 7-ray event candidates after image cuts; then dif- 
ferent positions in the FoV are used as background control regions. An alternative 
method (the so-called 'template model' [1]) uses hadron candidates which were 
reconstructed to originate from the same direction as the source candidate. In 
this normalisation is needed to estimate the background level in the 7-ray 

regime from the event counts in the hadron domain; this value must be derived 
from an average across the whole FoV. 

Both background estimates rely on the homogeneity of the (7-ray or ha- 
dron) acceptance across the FoV. While detector acceptance inhomogeneities are 
typically of the order of 5% or less, they may reach 10-20% in special cases such 
as large zenith angle observations or strong sky brightness variations due to stars. 
Much effort went into the development of an acceptance correction which accounts 
for as many as possible known systematic effects. Figure 1 shows as an example 
the acceptance at the positions of stars in the FoV relative to the surrounding 
acceptance, as a function of the B-band of stellar magnitude. The values were 
derived from data, averaging over all stars in the respective magnitude interval 
which have no other stars nearby. The lines show empirical parametrisations 
which are used as part of the overall acceptance correction. The functional de- 
pendence for the 7-ray and hadron regimes is quite similar, supporting the idea 
that the 'template' background estimate may be better able to cope with in- 
homogeneities (of whatever type) at individual pointings. Fortunately, the two 
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HEGRA lACT SYSTEiyi - EXPOSURE MAP 



Fig. 2. Exposure sky map of the HEGRA TACT system for the years 1997 to 2001, 
in Galactic coordinates. Darker colour denotes longer observations. 



different discussed background estimates are affected by different uncertainties of 
the acceptance determination; since the results derived with both methods are 
generally in good agreement, we conclude that both estimates are reliable. 

The 'template' background estimate is in principle applicable to the inves- 
tigation of source candidates which cover a large portion of the FoV, or to look 
for diffuse 7-ray emission. However, the cited systematic problems in the FoV ho- 
mogeneity prevent us currently from searching for source candidates which cover 
a large portion of the FoV. 

3. Search for new gamma-ray sources 

In the current analysis, two different grid search patterns were applied: a 
dense grid of 0.0625° spacing with a tight angular cut optimized for point sources, 
and a wider grid of 0.125° with a wider angular cut for slightly extended sources. 
In both cases, the significances for neighbouring grid points are correlated since 
the cuts are a bit larger than the grid spacing. Strong known 7-ray sources were, 
a priori, excluded from the search grid. 

Figure 3 shows significance distributions for the wider grid from a prelimi- 
nary analysis containing data up to 2001. The plot on the left hand side represents 
the search for new 7-ray sources. Here, 7-candidates in a ring surrounding the 
search position were used for background estimate. An excellent agreement with 
the Gaussian expectation for most of the scan positions is obtained, proving that 
the method works well. A few scan points show a 7-ray signal at a level of 5 — 7 u; 
these include objects that were detected as TeV emitters already in earlier analy- 
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flat-field check with hadron candidates 
signal search using gamma-candidates ( j q ^^^^^ higher statistics per search bin) 

Fig. 3. Significance distributions for the 0.125° grid. On the left panel, the result of 
the actual search for 7-ray sources is shown; in this case, 7-ray-like events in a ring 
surrounding the search position were used for the background estimate. The solid 
line indicates the expectation for a purely Gaussian distribution. On the right side, 
hadron-like events were treated in the same way as the 7-candidates to perform 
a check of the flat-field. Since the event statistics is 10 times higher than in the 
7-ray regime, the nearly Gaussian behaviour of this distribution shows that the 
background estimate is well under control. 



ses which were dedicated to these source candidates (e.g. CasA, H 1426+428). A 
few new 7-ray source candidates have been identified; final results for the entire 
data set are under way and will be presented at the conference. Strong unknown 
TeV emitters in the sky - at least where HEGRA was pointed - can however be 
excluded. 

The analysis was also performed for hadron events where no localized 
sources are expected. The hadron domain has ten times more counts per search 
bin than the 7-ray regime, the significance distribution (right panel of Fig. 2) 
is hence a very sensitive tracer for possible systematic errors in the fiat-field 
procedure. The distribution is nearly Gaussian, which is a good additional check 
that systematic acceptance variations are well under control. 
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Abstract 

Between the beginning of 1997 and fall 2002, the HEGRA collaboration 
was operating a stereoscopic system of 4 (later 5) Imaging Atmospheric Cherenkov 
Telescopes (lACT) on the Canary Island La Palma. We review the calibration 
schemes which were developed for the system, and report on the performance of 
the detector over this period. The system had an energy threshold of 500 GeV 
under optimum detector conditions at zenith. With the calibration schemes de- 
scribed here, a systematic accuracy of 15 percent on the absolute energy scale 
has been achieved. The continuous sensitivity monitoring provides a relative ac- 
curacy of a few percent, and shows that the threshold did not exceed 600 GeV 
throughout the whole period of operation. The readout electronics and the imag- 
ing quality of the dishes were well monitored and stable. The absolute pointing 
had an accuracy of better than 25". 

1. Introduction 

The HEGRA stereoscopic Imaging Atmospheric Cherenkov Telescope sys- 
tem was located in the Canary Islands, at 2200 m above sea level on the Roque 
de los Muchachos on La Palma (17°52'34" West, 28°45'34" North). The system 
consisted of 5 identical telescopes (CT2 - CT6), which operated in coincidence 
for the stereoscopic detection of air showers in the atmosphere. Now, after the 
disassembly of the system, we present an account of the long-term performance 
and stability of the system. A far more detailed description can be found in [2]. 

2. Telescope pointing and point spread function 

The 3.4 m reflectors were azimuthally mounted and driven by stepper mo- 
tors. The telescopes' pointing was calibrated and monitored offline by dedicated 
calibration runs, so called point runs, which used stars as reference objects. Those 
calibration runs were typically performed every few months. During data taking, 
the tracking algorithm relied solely on the position of the telescopes' axes, which 
were measured by 14 bit optical shaft encoders. The pointing was corrected off- 
line using an analytical model of the mechanical structure of the telescopes, the 
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Fig. 1. Reconstructed center of TeV emission from the Crab Nebula, and the Chandra 
X-ray image for reference. The white cross denotes the result for the entire HEGRA 
data set, the insets show the deviation from the expected position in yearly time 
intervals. 



parameters of which were determined from the point runs. The resulting pointing 
accuracy was 25" in right ascension and decHnation. 

The temporal stability and accuracy of the pointing calibration was verified 
for example by the examination of the center of gravity of the TeV emission of 
the Crab Nebula. In Figure 1, the white cross denotes the reconstructed center of 
the TeV emission for the whole HEGRA data sample; for reference, the Chandra 
X-ray image was placed at the year 2000 coordinates of the pulsar. Additionally, 
the HEGRA data were split up into yearly intervals, with results as shown in the 
insets. Within the systematic and statistical errors, the centers of the X-ray and 
TeV 7-ray emission are in agreement over the full lifetime of the experiment. 

Each telescope dish held 30 spherical glass mirrors with diameters of 60 cm; 
the total mirror area per telescope was 8.5 m^. The quality of the optical point 
spread function (psf) was monitored using the same point runs. The influence 
of the psf on 7-ray induced shower images is best obtained from the regular 
observations of strong TeV 7-ray sources. For two short periods only, the mirror 
alignment had been distorted by ice such that the image shape parameters needed 
to be tuned, until the mirror alignment was corrected. For the rest of the system 
lifetime, the image parameters were stable. 
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Fig. 2. HEGRA system sensitivity results, (a) cosmic ray event rate from runs below 
30° zenith angle (values rescaled to a 4-telescope setup); rates far below the expec- 
tation were taken under poor weather conditions, (b) electronic (kci) ^-nd relative 
optical (n'^pt) efficiency, averaged over all telescopes, (c) absolute optical efficiency, 
as obtained from muon runs, (d) Gamma-ray sensitivity of the HEGRA system 
(4-telescope setup), derived from Crab observations. 



3. Camera electronics calibration, relative sensitivity monitoring 

Each telescope held a camera with 271 photomultipliers, viewing 4.3° of the 
sky. The signals were digitized with 120 MHz FADCs. Less than 2% of channels 
were excluded due to technical problems in 95% of the runs. The relative gain 
of the camera pixels and the relative pixel timing was calibrated by laser runs, 
in which the camera was illuminated uniformly by laser flashes. These runs were 
performed every night, to provide a continuous FADC gain fiat-field. In addition, 
the individual pixel HV's were adjusted every few months, based on the laser run 
results; a trigger acceptance fiat-field of 5% or better was thereby achieved. 

Continuous monitoring of changes in the light sensitivity of the complete 
detectors was provided by the event trigger rate which originates from the steady 
background flux of charged cosmic rays (CR) (see Fig. 2a). The numbers directly 
translate into a change of the energy threshold of the system. Assuming an energy 
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threshold of 500 GeV under optimum conditions (see next section), the threshold 
deduced from the CR rate never exceeded 600 GeV. 

The main contribution to the sensitivity decrease of the detectors came 
from the electronics chain, presumably due to aging of the PM's last dynodes. 
The absolute gain of the camera electronics was monitored using the laser 
runs, values are given in Fig. 2b. To compensate for these losses, the global HVs 
of all cameras were increased three times (indicated by the vertical dashed lines) . 
The remaining part of the sensitivity loss was caused by losses in the optical 
throughput, presumably mirror aging (cf. k'^^^ in Fig. 2b, and also Fig. 2c). 

4. Absolute energy threshold, gamma-ray sensitivity 

The absolute energy threshold E'thr of the system was obtained by shower 
and detector simulations [1]; for early detector conditions, a value of Ethr = 
500 GeV ± 15%(sys) was estimated. The standard method of comparing the CR 
event rate to expectations from simulations is limited by uncertainties in the CR 
energy spectrum and composition. However, the detector sensitivity itself can 
be derived from the investigation of muon images in the cameras, and also from 
special laser calibration setups using a calibrated photodiode as reference. Figure 
2c shows results from muon runs, which were performed a few times per year; 
the values already include the correction by Hd.tei- The muon run results agree 
within 15% with the values used for the CR event rate simulations, confirming 
the estimated energy threshold. 

Finally, in Fig. 2d we show the sensitivity of the HEGRA system for weak 
7-ray sources, derived from Crab observations. The values denote the expected 
significances for a ten hour observation, in units of standard deviations a, for a 
7-ray source which has 10% of the Crab's source strength. The slight sensitivity 
decrease over the detector lifetime is presumably caused by a deterioration in the 
background rejection power, which could not be compensated for by the global 
HV increases [2]. 

To conclude, the telescope system was very well understood, regarding 
both the absolute calibration and the slight performance changes over the years. 
We believe that with the HEGRA system, an important contribution has been 
made to the effort of establishing ground based imaging Cherenkov telescopes 
as precision detectors, well suited within the broad range of astronomical instru- 
ments. 
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Abstract 

The first unidentified TeV source in Cygnus is confirmed by follow-up 
observations carried out in 2002 with the HEGRA stereoscopic system of Cerenkov 
Telescopes. Using all ~279 hrs of data, this new source TeV J2032-|-4130, is 
steady over the four years of data taking, is extended with radius 6.2', and has a 
hard spectrum with photon index —1.9. Its location places it at the edge of the 
core of the extremely dense OB association, Cygnus 0B2. Its integral flux above 
energies E > 1 TeV amounts to ~3% of the Crab nebula flux. No counterpart at 
radio, optical and X-ray energies is as-yet seen, leaving TeV J2032+4130 presently 
unidentifed. Summarised here are observational parameters of this source and 
brief astrophysical interpretation. 

1. Introduction &; Data Analysis 

Analysis of archival data (~ 121 h) of the HEGRA system of Imaging 
Atmospheric Cerenkov Telescopes (HEGRA lACT-System see for e.g. [13]) de- 
voted to the Cygnus region revealed the presence of a new TeV source [1]. This 
serendipitous discovery is now confirmed in follow-up observations from 2002 
(~158 h) by the same telescopes. Given the lack of a counterpart at other energies, 
TeV J2032-I-4130 may represent a new class of particle accelerator. 

For these analyses, cosmic-ray (CR) background events are rejected via a 
cut on the mean- s caled-width penameter w [2]. Event directions are reconstructed 
using the so-called 'algorithm 3' [6], and a cut is made on the angular separation 
6 between the reconstructed event and assumed source direction. So-called tight 
cuts are implemented: 9 < 0.12°, w < 1.1, and also demanding a minimum 
'^tei > 3 images for the 9 and w calculation. The background is estimated using the 
template model [1,14], and consistent results are also obtained using an alternative 
displaced background model which employs ring-segments within the FoV. For the 
centre of gravity (CoG) and source extension determination, an additonal cut on 
the estimated error in reconstructed direction (e < 0.12°) is applied, reducing 
systematic effects (e.g. [7]). The CoG and source extension are estimated by 
fitting a 2D Gaussian convolved with the instrument's point spread function to a 
histogram of 7-ray-like {w < 1.1) events binned over a 1° x 1° FoV. At the CoG, 
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Table 1. Numerical summary for TeV J2032+4130 (preliminary), (a) Centre of Grav- 
ity (CoG) and extension Usrc (std. dev. of a 2D Gaussian); (b) Event summary. 
The values s and b are event numbers for the 7-ray-like and bakground (from the 
Template and Displaced models, see text) respectively, and s — ab is the excess using 
a normalisation a. S denotes the excess significance using Eq. 17 of [9]; (c) Events 
after spectral cuts; (d) Fitted power law. 



(a) CoG & Extension (e < 0.12°) 



RA 02000 : 
Doc (52000 : 

O'src 



20'" 31™ 57.0= 
41° 29' 56.8" 
6.2' 



±1.2^,^, ±0.9^y. 



(b) Tight cuts: < 0.12°, w < 1.1, THel > 3 



Backer. 



S 



Template 
Displaced 



1245 
1245 



5926 
15492 



0.168 
0.065 



252 
243 



+7.1 
+7.1 



(c) Spectral Cuts: Tight Cuts + core< 200m 



Backgr. 



a b 



Energy estimation method: See [8] 

Displaced 974 5122 0.143 242 +7.9 



(d) Fitted Spectrum: Pure Power-Law 



dN/dE = Af (£;/lTeV)-T pli cm-^ g-i TeV" 

N = 5.3 (±2.2atat ± 1.3aya) X lO-^^ 

7 = 1.9 (±0.3stat ±0.3sys) 



the excess significance now exceeds 7a from all 278.2 h of data, and the source 
extension is confirmed as non point-like. The reconstruction of event energy 
follows the method of [8] and we use tight cuts plus a cut on the reconstructed 
air-shower core distance of the event core<200 m. A pure power law explains 
well the energy spectrum, showing no indication for a cut-off. A lower limit to 
the cut-off energy Ec ~ 3.6, 4.2 and 4.6 TeV is however estimated when fitting 
a power law+exponential cutoff term exp{—E/Ec) and fixing the power index 
at values 7 =1.7, 1.9 and 2.2 respectively. The integral flux for energies E > 1 
TeV is 5.9 (±3.1stat x 10"^'^), ph cm~^ s~^ or about 3% of the Crab nebula flux. 
Results are summarised in Table ITI and Fig [T] (upper panel). 

2. Modelling TeV J2032+4130 

Possible origins of TeV J2032+4130 have been discussed in literature [1,4,12]. 
One interpretation involves association with the stellar winds of member stars in 
Cygnus 0B2, individually or collectively, which provide conditions conducive to 
strong and stable shock formation for particle acceleration. Certainly the exis- 
tence of TeV emission suggests particles accelerated to multi-TeV energies. We 
have therefore matched the spectral energy distribution of TeV J2032+4130 with 
coarse leptonic and hadronic models (Fig. ITI lower panel). Another scenario in- 
volves particle acceleration at a termination shock, which are expected at the 
boundary where a relativistic jet meets the interstellar medium. TeV J2032+4130 
actually aligns well within the northern error cone of the bi-lobal jet of the famous 
microquasar Cygnus X-3 [10,11]. 

For simplicity we assume the TeV emission arises from either a pure sam- 
ple of non-thermal hadronic or leptonic parent particles. Under the hadronic 
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Cygnus OB2 Field: HEGRA CT-System 
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Fig. 1. Upper: Skymap of event excess significance (a) from all HEGRA CT-System 
data (3.0° x 3.0° FoV) centred on TeV J2032+4130. Some nearby objects are 
indicated (GeV sources with 95% contours). The TeV source centre of gravity 
(CoG) with statistical errors, and error circle (65% confidence) for the extension 
(std. dev. of a 2D Gaussian, agrc) are indicated by the white cross and white circle 
respectively. Lower: Spectrum of TeV J2032-I-4130 (labelled HEGRA) compared 
with purely hadronic (Protons E< 100 TeV) and leptonic (Electrons E< 40 TeV) 
models. Upper limits, contraining the synchrotron emission, are from the VLA and 
Chandra [3] and ASCA [1]. EGRET data points are from the 3rd EGRET catalogue 
[4]. 
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scenario the 7r°-decay prediction explains well the TeV flux when using a parent 
proton power law spectrum of index —2.0 with a sharp limit up to energies 100 
TeV. The neighbouring EGRET source 3EG J2033+4118 (likely not related to 
TeV J2032+4130) provides no constraint on this model. Associated synchrotron 
X-ray emission would also be expected from tertiary electrons (tt^ . . . /i^ . . . — >■ 
. . .). We have not yet modelled this component which essentially represents 
an absolute lower limit on any synchrotron emission visible. Assuming a pure 
leptonic scenario, TeV data are matched well by an inverse-Compton spectrum 
(up-scattering the cosmic microwave background) arising from an uncooled elec- 
tron spectrum with power law index ~ —2.0 and hard cutoff at 40 TeV. This 
allows us to predict the synchrotron emission as a function of local magnetic field, 
constrained by the available upper limits at radio and X-ray energies. The most 
conservative synchrotron prediction arises from the Bq = S.OfiG choice, which 
is realistically the lowest such field expected in the Galactic disk. But in fact, 
much higher fields {Bq > lO/iG) are generally expected in such regions containing 
young/massive stars with high mass losses and colliding winds (e.g. [5]). Deep 
observations by XMM and Chandra will provide strong contraints on the leptonic 
component. 
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Abstract 

Results from a search for TeV 7-ray emission from the vicinity of five X- 
ray-studied GeV sources is presented. A number of possible X-ray counterparts 
have been suggested based on ASCA and ROSAT observations (e.g. the SNR 
CTA 1, and X-ray binary LS I-|-61°303). Our search has yielded no convincing 
evidence for TeV emission at these GeV source positions and also a number possi- 
ble counterparts. Preliminary upper limits in the range ~1% to 10% of the Crab 
flux have been estimated above various energy thresholds 0.7 to 1.3 TeV. 

1. Introduction &: Data Analysis 

A number of EGRET sources visible above 1 GeV [8,10] were earlier stud- 
ied in X-rays with ASCA [12], revealing possible counterparts. The large field of 
view (FoV) of the HEGRA lACT-System (Imaging Atmospheric Cerenkov Tele- 
scope System, see [11]) allows good coverage of many GeV sources even if they 
were not the original targets of observation. We have selected five of these GeV 
sources for further analysis: 

GeV J0008-|-7304: This source could be associated with the nearby (~ 0.2° 
distant) supernova remnant CTA 1 (radio shell 100' diameter) and may con- 
tain a radio quiet pulsar [3]. This pulsar may power the extended, centre-filled 
non-thermal X-ray emission from CTA 1. Our observations were centred on the 
ROSAT point source RX J0007.0+7302 [15]. Previous TeV observations by the 
CAT [7] and Whipple [5] collaborations give upper limits for CTA 1. 
GeV J0241-I-6102: The GeV contours are quite consistent with the unusual X- 
ray binary LS I-|-61°303 [4]. An association with the COS B source 2CG135-I-01 
in the past has led to various 7-ray production models for LS I+61°303. Our 
observations were centred on 2CG135+01. 

GeV J1907-|-0557: This GeV source is located ~ 2° from the microquasar 
SS-433. ASCA studies reveal a weak source possibly comprised of two pointlike 
components within the GeV 95% contour. Our observations were taken from the 
HEGRA SS-433 archive. 

GeV J2026-|-4124: ASCA studies reveal a single point source within the 95% 
contour. This GeV source is ~ 1° from Cyg X-3. 
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GeV J2035+4214: Two pointlike (possibly non-thermal, labelled Srcl & Src2) 
and one marginally extended (thermal, Src3) ASCA sources are seen within the 
GeV 68% contour. HEGRA observations of these latter two GeV sources come 
from the extensive archive on Cyg X-3 and TeV J2032-I-4130 [14]. 

In these analyses we employ the mean- scaled-width w < 1.1 [2] cut as a 
means to reject the cosmic-ray background, and also a cut in the angular sep- 
aration 6 between reconstructed (employing the so-called 'algorithm 3' [6]) and 
assumed arrival directions for each event. A point-source search utilises a ^cut < 
0.12° cut. Extended sources with radius agrc utilise a cut 6'cut < \/^src + 0.12^ . A 
minimum of three telescope images per event ritei > 3 are demanded in calculating 
w and 6. The cosmic-ray background is estimated using the so-called template 
model [1,13]. Quite consistent results are also obtained using independent (e.g. 
so-called displaced) background models. 

Table 1. Numerical results for various GeV sources and other positions of interest. 
The values s and b are respectively the source counts and background estimate from 
the template model. 

Source Obs. time O^t s b ^ S'^ ^^99^ 

[h] [TeV] [deg] H 

GeV J0008+7304 



GeVJ0008+7304 


26.0 1.3 0.120 


72 


615 


0. 


114 


-1-0.2 


0.74 


RXJ0007.0+7302 


26.0 1.3 0.120 


70 


598 


0. 


114 


-1-0.2 


0.73 


CTA V 


26.0 1.3 0.195 


181 


1618 


0. 


114 


0.2 


1.09 




GeV J0241+6102 














GcVJ0241+6102 


28.8 0.8 0.120 


55 


338 


0. 


130 


+1.5 


1.74 


LS I+61°303 


28.8 0.8 0.120 


59 


380 


0. 


130 


+1.3 


1.50 




GeV J1907+0557 














GcVJ1907+0557 


114.1 0.7 0.120 


183 


1019 


0. 


164 


+1.1 


0.66 


AXJ1907.4+0549 


114.1 0.7 0.120 


187 


1101 


0. 


164 


+0.4 


0.62 




GeV J2026+4124 














GeVJ2026+4124 


275.0 0.7 0.120 


541 


3429 


0. 


168 


1.3 


0.29 


AXJ2027.6+4116 


275.0 0.7 0.120 


771 


4497 


0. 


168 


+0.6 


0.42 




GeV J2035+4214 














GcVJ2035+4214 


275.0 0.7 0.120 


784 


4457 


0. 


168 


+1.2 


0.45 


AXJ2036.0+4218(Srcl) 


275.0 0.7 0.120 


620 


4005 


0. 


168 


-1.8 


0.22 


AXJ2035.4+4222(Src2) 


275.0 0.7 0.120 


663 


4083 


0. 


168 


0.7 


0.27 


AXJ2035.9+4229(Src3) 


275.0 0.7 0.120 


584 


3626 


0. 


168 


0.9 


0.28 



a. Estimated 7-ray threshold energy for mean zenith z: E = 0.5cos(z)~^'^. 

b. Normalisation factor for the background b. 

c. Statistical significance using Eq. 17 of [9] . 

d. (jy'Jf^'' =99% upper limit E > _Eth [xlO^^^ ph cm-^s"!] 

e. For CTA 1, a radius cfsTc = 0.153° is used (the inner region for ROSAT spectral analysis [15]). 
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GeV J0008+7304 GeV J0241+6102 




R A (J2000) [ hrs ] R A (J2000) [ hrs ] 



GeV J2035+4214 




Fig. 1. Preliminary 2D skymaps of ex- 
cess significance S (only S > 0.5a are 
shown) for 2° x 2° FoV centred on 
the GeV sources (99% ellipse given). 
Events are integrated in a circle of ra- 
dius Ocut < 0.12° centered at positions 
over an 0.025° x 0.025° grid. Other po- 
sitions of interest according to the text 
are shown. In the GeV J2026-h4124 and 
GeV J2035-t-4214 skymaps, the uniden- 
tified TeV source TeV J2032-F4130 is 
clearly seen, somewhat 'clipped' since a 
maximum S = +Aa is set. 
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2. Results &: Discussion 

As a-priori targets, we have chosen the ASCA sources identified by [12], 
some pre-identified objects (e.g. CTA 1, LS I+61°303 etc.), and also a point 
source centred on each GeV source (SIMBAD coordinates). Numerical results are 
presented in Table ITT] and 2D skymaps of excess significances are given in FiglTlFor 
all of our a-priori-chosen positions, no convincing evidence for TeV emission was 
found. Upper limits at the 99% confidence level were estimated for each position 
using the expected Crab event rate from Monte Carlo simulations scaled according 
to an estimate of the 7-ray acceptance in the FoV. These preliminary upper limits 
are in the range 1% to 10% of the Crab fiux at various energy thresholds (i?th ~ 

0. 7 to 1.3 TeV, based on the mean zenith angle of events). Further discussion and 
comparison with models will be presented at the conference. 
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Abstract 

In total 54 Active Galactic Nuclei (AGN) apart from Mrk-421 and Mrk- 
501 have been observed between 1996 and 2002 with the HEGRA Cherenkov 
Telescopes in the TeV energy regime. Among the observed 54 AGN are the 
meanwhile well established BL Lac type objects H 1426+428 and lES 1959+650. 
The BL Lac object lES 2344+514 and the radio galaxy M 87 show evidence for a 
signal on a 4(T level. Fluxes resp. upper limits are given for each of the 54 AGN. 

1. Introduction 

In the commonly adopted view AGN are powered by a central super mas- 
sive black hole with ~ 10^ Mq surrounded by an accretion disk. Two relativistic 
plasma outflows (jets) perpendicular to the accretion disk are pointing in oppo- 
site directions away from the center black hole [1]. Most detections of extra- 
galactic TeV 7-rays refer to objects of the BL Lac type, i.e. AGN having their 
jet pointing close to the observers line of sight. In this work the results of an 
analysis of the observations of 54 AGN with the HEGRA telescope system are 
presented. The HEGRA collaboration was operating a stereoscopic system of 
imaging air Cherenkov telescopes (lACT-System) [2] and one stand alone tele- 
scope [3] (not used for the present analysis) on the Canary island of La Palma 
(28.75° N, 17.90° W) at an altitude of 2200 m a.s.l. The stereoscopic reconstruc- 
tion technique first introduced by HEGRA allows the complete reconstruction of 
the shower geometry and results in improved angular and energy resolutions as 
well as a very good 7-hadron separation (hadron rejection up to a factor of 100) 
using the so called mean scaled width (mscw) parameter. The data set used for 
this work amounts to a total exposure time before rejection of runs with poor 
quality of ~ 1150 hours (without Mkn-421 and Mkn-501) corresponding to more 
than one year of continuous observations in moonless nights with the HEGRA 
lACT-System. After quality cuts, 1017 hours remain in the analysis. The applied 
event cuts are optimized on data of the Crab Nebula which is used as a calibration 
source in TeV astronomy. The significance of an excess is calculated following [4] . 
Observed Crab-7-rates are used to compute integral flux upper limits [5]. 
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, I I- 17.51 / 12 
Constant 5.311 ± 0.9708 
Mean -0.9879E-01 ± 0.1560 
Sigma 1.031 ± 0.1261 

□ distribution of AGN significances 
Gaussian fit (-2.5o to +2.5o) 



1 

"5 5 10 15 20 25 

DC significance S, [a] 

Fig. 1. Distribution of significances for all AGN analyzed in this work. The objects 
lES 1959+650 (1) and H 1426+428 (2) show a clear deviation from the background 
expectation, represented by a Gaussian fit from -2.5 a to 2.5 a. Further two objects, 
lES 2344+514 (4.4 cr, (3)) and M87 (4.1 a, (4)) also show a deviation from the 
background expectation and thus for emission of TeV 7-radiation. 

2. Results 

A distribution of the significances for steady state emission (DC) of all 
analyzed objects is shown in Fig. 1. The objects lES 1959+650, H 1426+428, 
lES 2344+514 and M87 show significances deviating from the expectation of a 
Gaussian distribution in case of pure background fluctuations. In Tab. 1 a list 
of all objects analyzed in this work ordered by ascending redshift is shown along 
with their observation time and the results of this analysis. The results of the 
observations of the BL Lac objects lES 1959+650 and H 1426+428 as well as a 
tentative detection of a signature of absorption by pair production on the diffuse 
extragalactic background radiation in the TeV spectrum of H 1426+428 were re- 
ported elsewhere [6,7]. The results of a dedicated analysis of the M87 data are 
presented separately at this conference and have been published recently [8,9]. 
The first TeV detection of lES 2344+514 was reported by the Whipple group in 
1998 with 6(7 in one night only [10]. The results of the HEGRA observations of 
lES 2344+514 of the year 1997 and 1998 were first reported by [11] with a DC 
significance of 3.3 ex. The analysis presented here includes the complete dataset 
(1997 - 2002) and results in an excess of 63 ± 14 7-like excess events (Non = 235, 
(Noff) = 171) on a significance level of 4.4 a and no evidence for burst-like be- 
haviour. However, the excess is accumulated almost exclusively in the observation 
periods of August to November 1998 and September 2002 (61 ±12 excess events). 
In Fig. 2. the distribution of the squared angular distance of the reconstructed 
shower events to the object position as well as a plot of the excess rate vs. obser- 
vation periods are shown. The observed excess results in a flux of <I'(E>0.80 TeV) 
= (0.8 ± 0.3) ■ 10~^^ cm~^ s~^ corresponding to 3.3% of the Crab-Nebula flux. 
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Fig. 2. Left: Distribution of reconstructed squared angular distances of the 
lES 2344+514 data. The distribution of the on-source events is represented by 
the data points. The background is shown as a shaded histogram. Right: The 
excess rate in 7/hour for the three observation campaigns on lES 2344+514. The 
excess is accumulated almost exclusively in 1998 and 2002. 



3. Summary 

The two objects lES 1959+650 and H 1426+428 can be considered as well 
established sources of TeV 7-radiation. The evidence for lES 2344+514 being 
a TeV source is strong taking the detection by the Whipple collaboration into 
account. The detection of TeV 7-rays from M87, if confirmed, would be the first 
detection of photons in the TeV energy regime of an AGN not classified as a BL 
Lac object. 
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Table 1. Results of all objects of the HEGRA AGN data sample. The object names 
are given as well as observation time, number of ON- and OFF-events, energy 
threshold, upper limits in Crab units and in flux units and fluxes in Crab units. 



Object 


z 


time 


Non 


Noff 


Efhr 


FTTT^'^f> B+v.r) 


$^7°(> Et,„) 
[lO-^^^cm-^s-i] 


F(> Ethr) 






[h] 


# 


# 


[TeV] 


[Crab] 


\Crah\ 


lES 0647+250 





4.1 


12 


10 


0.74 


0.13 


3.65 




MG 0509+0541 




15.8 


60 


50 


0.79 


0.11 


2.61 




M87 


.004 


70.0 


241 


184 


0.76 






0.033 


NGC 315 


.016 


14.6 


31 


33 


0.75 


0.05 


1.29 




NGC 1275 


.018 


87.6 


231 


236 


0.75 


0.03 


0.84 




H 1722+119 


.018 


5.1 


21 


14 


0.76 


0.21 


5.51 




PKS 2201+04 


.028 


17.8 


59 


46 


0.79 


0.08 


1.89 




V Zw 331 


.029 


4.1 


9 


9 


0.75 


0.09 


2.44 




NGC1054 


.032 


57.9 


134 


155 


0.75 


0.02 


0.46 




3C 120 


.033 


25.4 


64 


70 


0.78 


0.05 


1.14 




NGC 4151 


.033 


7.0 


16 


18 


0.74 


0.07 


1.99 




UGC01651 


.037 


14.3 


44 


36 


0.74 


0.07 


1.79 




UGC03927 


.041 


6.3 


7 


16 


0.86 


0.09 


1.93 




lES 2344+514 


.044 


72.5 


235 


171 


0.80 






0.033 


MknOlSO 


.046 


9.8 


29 


32 


1.38 


0.12 


1.24 




lES 1959+650 


.047 


163.7 


1202 


454 


1.17 






0.053-2.200 


3C 371.0 


.050 


5.4 


16 


18 


1.40 


0.19 


1.88 




B2 0402+37 


.054 


6.7 


9 


17 


0.74 


0.05 


1.32 




I Zw 187 


.055 


16.0 


44 


32 


0.78 


0.09 


2.21 




Cyg-A (3C 405.0) 


.057 


59.0 


159 


161 


0.77 


0.03 


0.83 




IBS 2321+419 


.059 


22.3 


53 


66 


0.76 


0.03 


0.86 




30 192.0 


.060 


2.9 


8 


7 


0.78 


0.20 


5.01 




40+31.04 


.060 


3.0 


8 


9 


0.74 


0.14 


4.01 




BL Lacertae 


.069 


26.7 


94 


67 


0.87 


0.28 


5.97 




lES 1741+196 


.083 


10.2 


28 


26 


0.78 


0.07 


1.88 




40+01.13 


.084 


7.7 


30 


31 


0.81 


0.10 


2.45 




PKS 2155-304 


.116 


1.8 


4 


4 


4.62 


0.27 


0.39 




lES 1118+424 


.124 


2.0 


5 


4 


0.79 


0.24 


5.88 




IBS 0145+13.8 


.125 


3.2 


2 


1 


0.75 


0.06 


1.73 




EXO0706. 1+5913 


.125 


33.7 


81 


85 


0.86 


0.06 


1.32 




H 1426+428 


.129 


258.5 


796 


624 


0.77 






0.024-0.060 


30197.1 


.130 


15.0 


22 


24 


0.79 


0.05 


1.17 




lES 1212+078 


.130 


2.4 


6 


8 


0.78 


0.17 


4.26 




IBS 0806+524 


.138 


1.0 


2 


2 


0.86 


0.29 


6.19 




IBS 0229+200 


.139 


3.0 


11 


8 


0.78 


0.17 


4.25 




RBS 0958 


.139 


3.8 


18 


9 


0.75 


0.28 


7.57 




IBS 1255+244 


.140 


5.9 


14 


14 


0.79 


0.12 


2.88 




MS1019. 0+5139 


.141 


17.5 


44 


43 


0.78 


0.07 


1.78 




IBS 0323+022 


.147 


14.3 


24 


32 


0.81 


0.04 


1.00 




OQ 530 


.152 


9.4 


32 


30 


0.89 


0.10 


2.05 




30 273.0 


.158 


12.2 


44 


46 


0.91 


0.09 


1.82 




IBS 1440+122 


.162 


13.1 


35 


41 


0.77 


0.08 


1.95 




PKS 0829+046 


.180 


18.0 


59 


55 


0.81 


0.06 


1.35 




PG 1218+304 


.182 


3.9 


7 


8 


0.75 


0.12 


3.22 




IBS 0347-121 


.185 


1.9 


13 


7 


1.34 


0.56 


5.90 




IBS 0927+500 


.186 


13.3 


30 


29 


0.79 


0.06 


1.44 




PKS 2254+074 


.190 


16.3 


44 


48 


0.77 


0.05 


1.27 




MS0317.0+1834 


.190 


2.7 


5 


6 


0.74 


0.12 


3.30 




IBS 1011+496 


.200 


2.0 


1 


3 


0.81 


0.11 


2.58 




IBS 0120+340 


.272 


18.9 


36 


44 


0.75 


0.04 


1.03 




2B 0414+0057 


.287 


4.5 


18 


15 


0.81 


0.13 


3.07 




85 0716+714 


.300 


1.7 


6 


4 


1.46 


0.38 


3.54 




30 066A 


.444 


1.3 


3 


3 


0.75 


0.17 


4.70 




PKS 0219-164 


.698 


1.7 


5 


10 


1.67 


0.27 


2.04 





HEGRA Contributions to the 28th International Cosmic Ray Conference 



49 



Study of the VHE Gamma Ray Emission from the AGN 
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Ring 6, Munich 80805, Germany 

Abstract 

The BL Lac object 1ES1959+650 has been observed by the HEGRA Col- 
laboration in the years 2000, 2001 and 2002. Here we report on the results ob- 
tained in 2002 with the standalone CTl telescope. While the source has only 
shown weak TeV activity until mid-May 2002, significant TeV emission has been 
observed later in 2002. During 2002 the source has been monitored for about 200 
hours, of which 30 hours data were taken under moonlight conditions. The signal 
from the 'dark' time 2002 data set showed a significance of 11a. Preliminary 
results on the light curve and the spectrum of the 2002 'dark' time data up to 
45° zenith angle are presented. 

1. Introduction 

The galaxy 1ES1959-I-650 has a strong emitting jet from the nucleus, ori- 
ented along the line of sight. It was predicted by [8] to be a TeV 7-ray emitter 
candidate from scaling its X-ray observed spectrum, and from studies based on 
the synchrotron-self Compton model [2]. It was detected for the first time as a 
TeV 7-ray emitter by the Telescope Array experiment [7], that reported a signal 
with a significance of 3.9(T after 56.7 hours of observation in 1998. 

In 2000 and 2001 the HEGRA collaboration observed 1ES1959+650 with 
the CT System for 94 hours. The source was at low state and an averaged 
signal of 5.4o" of significance level was recorded [3]. In May 2002 the VERITAS 
Collaboration reported a strong TeV 7-ray signal from 1ES1959+650 [11]. The 
detection had a significance of 13(j from observation in 2 consecutive nights. The 
followed observations by the HEGRA Collaboration confirmed 1ES1959+650 as 
a highly variable VHE 7-ray source [1,4]. The results of observations with the 
HEGRA CT System are reported elsewhere on this conference. 

2. Observations with the HEGRA CTl Telescope and Data Analysis 

The HEGRA Cherenkov Telescopes are installed at the Roque de los 
Muchachos Observatory on the Canary Island of La Palma (28.75° N, 17.89° 
W, 2200 a.s.L). The telescope CTl [6] has a mirror area of 10.3 m^ and a camera 
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Table 1. Observation times of 1ES1959+650 with HEGRA CTl 



Observation conditions 



2000 2001 



2002 



hvOO 
hv04 
hv08 
hvl2 



nom 



good 'dark' night conditions, nominal HV 
weak moonhght, no HV reduction 



58 h 



86 h 
24 h 



156 h 
8.1 h 
16.1 h 
6.1 h 
2.9 h 



moonhght, 4% HV reduction 
moonhght, 8% HV reduction 
moonlight, 12% HV reduction 



12 h 



8 h 
7 h 



of 127 pixels of 0.25° diameter each. The camera has a FOV of 3° diameter. In 
zenith position, the telescope has an energy threshold of 750 GeV and a typi- 
cal source discovery sensitivity of 3.5 a X\/t, with t in hours, for a Crab-like 
source. Table 1 summarizes the CTl observation times of different background 
light conditions of 1ES1959+650. The observation angle ranged between 36.4° 
and 54.2° (with rather few data above 45°). A total of 1.6 million triggers have 
been recorded in 2002. 

Quality cuts are applied to the raw data in order to reject accidental noise 
triggers and nights with bad atmospheric conditions. The image parameters of 
the recorded events were calculated. The adopted 7/hadron separation method 
is described in [5]: the so-called dynamical super-cuts in the parameters DIST, 
WIDTH and LENGTH are applied, in order to take into account variations of the 
image parameters with the zenith angle, the energy of the primary particle and 
the impact parameter of the shower (estimated mainly from the image parameter 
DIST). The ALPHA parameter distribution after the dynamical super-cuts for 
'dark' nights data, taken up to a zenith angle of 45°, is shown in Figure 1. The 
number of recorded 7 events has been estimated by fitting this distribution with 
a Gaussian function in the signal region 0° < ALPHA < 13° and a smooth 
background (BG) has been estimated by fitting the same data sample with a 
polynomial of second order (without linear term) in the range 20° < ALPHA < 
80°, where no signal is expected. The BG was then extrapolated down to ALPHA 
= 0°. In Figure lb the BG parametrization is verified for a pure BG sample. 

3. Results 

The CTl data taken in 2000 and 2001 (not shown here) reveal that the 
source was in a low state. The averaged estimated 7-ray flux above ITeV was 
between 5% and 6% of the Crab flux at the same energy range [9,10]. 

The ALPHA distribution analysis for the data taken during the year 2002 
for ~ 150 h of observation (Figure 1) shows a 7-ray signal from 1ES1959+650 
with a significance of ll.lcr, with 664 ± 40 excess and 1672 ± 37 BG events for 
ALPHA < 13°. The error takes the BG extrapolation uncertainty into account. 
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ALPHA plot 1ES1959+650 
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Fig. 1. Left: Distribution of the image pa- 
rameter ALPHA for the 2002 data; max- 
immn zenith angle: 45°. Right: AL- 
PHA distribution of a normalized sample 
of off-data, plotted here for comparison. 



1ES1 959+650 Differential Spectrum, CT1 2002 



Energy (TeV) 

Fig. 2. Differential spectrum of the 
AGN 1ES1959+650 (HEGRA CTl 
2002 data set, no moonlight, 
max.ZA 45°). The hue that fits the 
data is a power law with a fixed cut- 
off energy of 2.4 TeV. 



The TeV light curve of the year 2002 is shown in Figure 3 together with 
the 2-12 keV activity monitored by RXTE ASM [12]. Two periods of stronger 
activity were recorded in that year: one in May (MJD 52414), when the Crab flux 
was reached, and a flare in July (MJD 52460 and 52470). The first comparisons 
between the X-ray fluxes and the TeV fluxes show no evident correlation between 
the two considered energy ranges. The differential spectrum of the source is 
shown in Figure 2. It can be fitted with a power law function, with a spectral 
index a = 3.49 ± 0.19 {x^/d.o.f. = 9.1 /6). On the other hand, the Extragalactic 
Background Light should already affect the spectrum above a few TeV. We tried 
also a power law approach modified by an exponential energy cutoff. Using a 
cutoff of Ec = 2.4 TeV [3] resulted in a somewhat better fit with a = 2.54 ± 0.21, 
X^/d.o.f. =7.1 /6. The data quality is too poor above a few TeV to draw decisive 
conclusions. 

4. Discussion and Conclusions 

The AGN 1ES1959+650 has been observed with the HEGRA CTl tele- 
scope. In May and July 2002, a short period of high emissions in the TeV range 
was observed, followed by a gradual decrease. The analysis of the 150 'dark' hours 
of observations with good conditions show a signal with a significance of ll.lcr. 
The averaged spectrum can be described by an unbroken power law, but a cutoff 
around 2.4 TeV is equally likely. A simple attempt to study X-ray- TeV correlation 
showed no significant correlation. In particular, we note that the X-ray flux did 
not vary significantly around the TeV flare at MJD 52460-52470. It is expected 
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1 Light Curve 1ES1959, 2002. HEGRA CT1 
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Fig. 3. The light curve of the data taken with HEGRA CTl in 2002 (above), in 
absence of moonhght. Below: data taken with RXTE ASM in the same period. 

that next generation low threshold lACTs might allow detailed clarifying studies. 
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Observation of VHE Gamma Rays from the Remnant of 
SN 1006 with HEGRA CTl 
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1. Abstract 

The HEGRA collaboration observed the shell-type remnant SN1006 with 
the telescope CTl for 219 h. For the observations we pointed the telescope towards 
the SNR center. An excess is observed from NE part of the rim with a statistical 
significance of ~ 5 a. Only large zenith angle (ZA) observations ( >70°) are 
possible from the telescope site. The energy threshold is ~ 18 TeV. We discuss 
the analysis and the problems related to large ZA observations. 

2. Introduction 

Supernova remnants (SNRs) are considered as the most likely sources of 
galactic cosmic ray. The shell type remnant G327.6+14.6 is originated from the 
SN 1006. This remnant is nearly circular (angular diameter 0.5°, radius of 
9.5 pc, distance 2.18 ± 0.08 Kpc )[1], limb-brightened and shows the NE and 
SW regions significantly brighter than the rest of the shell (bipolar or barrel- 
shaped) [1]. The shell expands with velocity 2890 ± 100 km/s [2] into a low density 
ISM. ASCA and ROSAT detected non-thermal X rays from the NE (the brighter 
one) and SW regions [5], [6]; while Cangaroo [7] detected TeV 7 rays only from the 
NE part of the rim. Non-thermal X rays were explained as synchrotron emission 
from 10-100 TeV electrons . The possible origins of the TeV radiation are : 1) the 
above electron population, which can produce TeV 7 rays via inverse Compton 
scattering on the CMBR and local IR photons ; 2) a possible contribution from 
accelerated hadrons (ttq => 77) [6] . 

3. Observations and Analysis 

We observed SN 1006 for 219 hours from 1999 to 2001 with the HEGRA 
telescope CTl. CTl is part of the HEGRA cosmic ray detector complex on La 
Palma (28.75° N, 17.9° W, 2225 m asl). SN 1006 (Ra 15:02:48.8, Dec -41:54:42) 
culminates at less than 20° elevation when observed from La Palma. Therefore 
only large zenith angle (ZA) observations are possible. We collected in total 
346000 so-called ON-source events restricting observations to a narrow zenith 
range between 71-73° ZA as well as a small set of OFF-source data at the same 
ZA. We recorded also a large background sample of muons passing through the 
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camera. Such muons leave narrow images, which might mimic 7 events. The 
data reconstruction proceeded in the following way: At first we applied the usual 
filter cuts rejecting accidental noise triggers or event samples recorded during 
non-optimal atmospheric conditions [8,9,10]. As a next step the image parame- 
ters were reconstructed and various cuts applied to select 7 events. Because of 
the large ZA observations we had to modify considerably our standard 7/h sep- 
aration procedure. The image dimensions shrink with the ZA (roughly oc ^^^l^zA) 
). The shape of small images is determined with less precision than in the zenith 
position because the pixel diameter is constant. Therefore the 7/h separation 
power decreases. The change of the image parameters could be either evaluated 
by Monte Carlo simulations or by studying real hadron events and assuming that 
7 image parameters would scale in the ratio. We opted for the second procedure 
and studied a large sample of hadronic images covering the entire ZA range be- 
tween 0° and more than 70° [10]. Up to 45° we could verify the scaling of the 
7 image parameters with events from a large sample of Crab Nebula and Mkn 
421 data. We assumed that the scaling is also applicable to data at 72° ZA. ( 
An error in the extrapolation should not fake a possible source but would mainly 
result in a less than optimal selection of 7s) . The scaled image cuts in Width and 
Length were applied to the data. Muons are more likely to generate background 
at large ZA observations because muons are more likely to pass simultaneously 
the glass envelope of many PMs. We studied the muon rate with closed camera 
and found a raw rate of 5.1x10^'^ Hz at 72° ZA. Efficient muon rejection cuts 
were developed [10] from comparing the measured muon event topology with MC 
simulated 7 candidates. Only a fraction of 1.8x10"'^ survived the cuts while the 
majority of 7s passed these cuts. 

For the search of a possible excess of localized 7 emission around SN 1006 
we applied the False Source Method (FSM). A grid of 0.1° step size was defined in 
the FOV and at each knot we searched for a pointlike 7 source. This method works 
also for slightly extended sources underestimating somewhat the 7 fiux. At each 
knot we applied the usual Distance cut and applied a 12.5° cut in the ALPHA 
distribution. The remaining sample is basically a coarse measurement of the 
hadron fiux across the FOV (it should be noted that there is a twofold ambiguity 
because one is not able to distinguish the head/tail of the shower). Fig lA shows 
the reconstructed hadronic background as well as a possible excess. Already at 
this step an excess around RA=15:04:10 and DEC =-41:40:00 is visible. In the 
next step we tried to subtract the background by extrapolating in each ALPHA 
distribution the background in the 0-12.5° ALPHA range. For this procedure we 
fitted a second degree polynomial (without the linear term) to the data between 20 
and 80°. Then we subtracted the extrapolated background into the low ALPHA 
regime from the actual data and determined the excess. The procedure has been 
extensively tested [9]. The excess map is shown in Fig IB. A clear excess is found 
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at the NE rim . Significance is in part influenced by the background statistics. 
We improved the background precision by averaging over a ring around the center 
of the FOV but leaving out a section, around the position under test , which 
might contain 7 events [10]. This method increases the background statistics by 
a factor ^ 3 and reflects the true behavior at the low ALPHA range. The lower 
panel in Fig 1 shows the ALPHA distribution of the excess bin together with the 
higher statistics background distribution in the ring. An excess of 103 ±17 over a 
background of 225 events has been determined using the recommended procedure 
[7]. The signiflcance of the excess is 5.6a . The excess location is nearly coincident 
with the observation of the Cangaroo group. If we allow for an uncertainty 9 bins 
we obtain a reduced signiflcance of 5. la. We verifled the correct positioning of 
the telescope at large ZA by measuring the position of two bright stars in the 
FOV. The positions of k Cen (mag=3.13) and (3 Lup (mag=2.68) allowed us to 
determine the excess position with 0.1° precision. For a 5a excess it is not possible 
to determine a spectrum. From MC simulations including all above mentioned 
cuts we calculated a threshold of 18 TeV and the collection area assuming an 
unbroken powerlaw with spectral index of -2.3 [5]. Using the experimental data 
and the parameters from the simulation we calculated a preliminary integral flux 
of $(>18 TeV) = (2.5± 0.5,tat)xl0-^3^. The evaluation for other photon 
indices and of the systematic error (at least 35%) is still ongoing. 

4. Summary and Conclusions 

From the analysis of the large zenith angle observation of SN 1006 we 
conclude: a) that we see evidence (5.1 a) for 7-emission at the NE rim of the 
shell coinciding within 0.1° with the Cangaroo observation at lower energies, b) 
The threshold of CTl at 72° ZA is 18 TeV assuming an unbroken powerlaw 
with a spectral coefficient of -2.3 . c) A preliminary integral flux of $(E>18 
TeV) = (2.5 ± 0.5sta4)xlO~^^^j^^ has been determined. The systematic error 
is still under study but is at least 35%. d) We observe no excess at the SW 
rim and determined an upper limit of 1.6x 10-^3 (^q.L. 99.5%) We conclude 
further that it is possible to carry out searches for 7 sources at large ZA although 
limitations from the pixel size and the large threshold are obvious. 
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Fig. 1. Fig lA is the sky map of ON-source data after 7 selection but before of sub- 
traction of residual background; Fig IB is the sky map of the excess (= ON-source 
- residual background); Fig IC is the sky map of excess significance; Fig ID is the 
distribution of the excess. Fig 1 lower panel is alpha distribution for ON-source 
data and background (ring method) at the NE rim position. All maps show the 
FOV of CTl; RA increases from right to left. The black circle in IB and IC shows 
the position of the SNR shell. The NE part of the rim is around the upper-left part 
of this circle. 



